Positronic complexes with unnatural parity 
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The structure of the unnatural parity states of PsH, LiPs, NaPs and KPs are investigated with 
the configuration interaction and stochastic variational methods. The binding energies (in hartree) 
are found to be 8.17 x 10~ 4 , 4.42 x 10~ 4 , 15.14 x 10~ 4 and 21.80 x 10" 4 respectively. These states are 
constructed by first coupling the two electrons into a configuration which is predominantly 3 P C , and 
then adding a p-wave positron. All the active particles are in states in which the relative angular 
momentum between any pair of particles is at least L = 1. The LiPs state is Borromean since 
there are no 3-body bound subsystems (of the correct symmetry) of the (Li + , e~ , e~ , e + ) particles 
that make up the system. The dominant decay mode of these states will be radiative decay into a 
configuration that autoionizes or undergoes positron annihilation. 

PACS numbers: 36.10.-k, 36.10.Dr, 34.85,+x 
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I. INTRODUCTION 

The stability of a bound state composed of two elec- 
trons and a positron, the positronium negative ion, was 
first demonstrated in a seminal calculation by Wheeler 
[l|. Shortly after this calculation, the four body systems, 
PsH and Ps2 were shown to be stable @, H| . Since that 
time, only a few other electronically stable states have 
been discovered that can be formed from combinations 
These are additional bound states of 



of P + ^&~ an d e 4 

Ps2 0, [E @, 0]) a compound that is best described as 
e+PsH @, and a (p+, 4e~, 2e+) complex 0. Addition- 
ally, a number of atoms have been identified as b eing 
capable of binding positronium and positrons 0, [ToL I 111 ] 
Just recently, a new class of positronic compounds that 
are electronically stable was identified |12|. The new PsH 
and NaPs bound states were unnatural parity states with 
symmetry conditions that act to prevent decay into the 
lowest energy dissociation products. An unnatural parity 
state is a state with parity equal to n = (—1) where 
L is the orbital angular momentum of the state. These 
PsH and NaPs systems have the two valence electrons in 
a spin-triplet state, a total orbital angular momentum of 
zero, and an odd parity, i.e. L n = 0~ . In addition, these 
states had the unusual feature of decaying very slowly by 
27 or 37 annihilation. 

In this paper, more details about the L n = 0~ nega- 
tive parity states of PsH and NaPs are given. Negative 
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parity states of LiPs and KPs are also identified as being 
electronically stable. The LiPs state has the additional 
distinction of being a Borromean state [H, [3, [H, Ell 
since the (Li + , e~, e~, e + ) system has no stable 3-body 
state that can act as a parent for the four body 2 ' 4 S° 
state. We have adopted the definition of Richard: A 
bound state is Borromean if there is no path to build it 
via a series of stable states by adding the constituents one 
by one (l3| . 

It should be noted that there are analogs of these states 
in the alkaline-earth sequence. Configuration interaction 
(CI) approaches has been used to demonstrate the sta- 
bilityof the Be", Mg", Ca" and Sr" np 3 4 S° states 
LTM| - However, the issue of whether an electron can 
be attached to the 3 P C state of H _ into an 4 S° state of 
H 2 ~ has been the subject of some controversy. A com- 
plex rotation method was applied to a large basis CI wave 
function and a shape threshold lying about 1.4 eV above 
the 3 P e threshold was predicted However, this was 
contradicted by a much more sophisticated hyperspheri- 
cal calculation that exhibited no sign of a resonance [2(| • 



II. THEORETICAL OVERVIEW 
A. Symmetry conditions for binding 

The stability of these systems lies in the symmetry re- 
lations between the pairs of particles that make up the 
system. The discussion of these conditions will be ad- 
dressed specifically to PsH, but these conditions, with 
some small modifications, will also apply to the other 
systems addressed in this paper. 

The electronic stability of PsH can be motivated by 
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consideration of the H~(2p 2 3 P C ) bound state [H|,[2|,[ll|. 
This state has an energy of —0.12535545 hartree [23j and 
is electronically stable due to symmetry conditions. It 
cannot decay into the H(ls) + e~ channel since the £ = 1 
partial wave of the electron automatically results in a 
state of negative parity. The L u = 0~ state of PsH is 
formed when the positron is trapped into a 2p state of the 
H~ attractive potential well. The possible decay modes 
are constrained by the symmetry conditions. Dissocia- 
tion into Ps(ls)+H(ls) is forbidden since II = (— 1) L 
where L is the orbital angular momentum between the 
Ps(ls) and H(ls) fragments. Similarly, dissociation into 
Ps(ns)+H(n£) or Ps(n^)+H(ns) does not occur since it 
is not possible to construct an L u = 0~ state if one of 
the angular momentum is zero. The lowest energy dis- 
sociation channel would be into Ps(2p)+H(2p) (with the 
two fragments in a p-wave) with an energy threshold of 
—0.1875 hartree. Another possible decay would be into 
the H~(2p 2 3 P e ) + e + channel but the threshold energy 
here is -0.12535545 hartree 

It is easy to see that there is potentially a large energy 
advantage associated with binding the positron to the 
negative ion. If the H~ state is regarded as a point par- 
ticle with an internal energy of « —0.125 hartree, then 
a positron in the 2p state will lower the total energy to 
—0.250 hartree. In actuality the H~(2p 2 3 P C ) state is 
diffuse 22], but the advantage of attaching the positron 
to the negative ion is clear. 



been calculated to be approximately 10 4 s _1 [13, H^|- 
Similarly, the 37 process, which happens when the anni- 
hilating pair are in a spin-triplet state, can also occur at 
a rate proportional to a 6 when the pair have a relative 
angular momentum of 1. The discussions below about 
the symmetry conditions for positron annihilation con- 
cern the fast 27 and 37 processes for pairs in relative 
s-states. 

Consider the electron-positron annihilation of a PsH 
state of 2 S° symmetry. The relative angular momentum 
of the annihilating pair (X re i) must be zero. This means 
the total angular momentum of the state will come from 
the center-of-mass motion of the annihilating pair (L cm ), 
and from the angular momentum of the spectator elec- 
tron (-^spectator)- The total parity of the state is deter- 
mined by the parity of the individual constituents, i.e. 
II = (-i) L =pc=tato,+Lcm+Lrc^ ft [ s simply not possible to 
form an odd parity state with a total angular momen- 
tum of zero if any one of the angular momenta is zero. 
Consequently, a two electron/one positron state of 2 S° 
symmetry cannot decay by the fast 27 process. 

These arguments also apply to the 37 annihilation pro- 
cess. The 37 process occurs for electron-positron pairs in 
a spin-triplet state with a relative angular momentum of 
zero. Once again, it is simply impossible to form a state 
of 2 S° (or 4 S°) symmetry if the relative angular momen- 
tum of the annihilating pair is zero. So it is reasonable 
to conclude that the lowest order 37 decay is not possible 
from a 2 ' 4 S° state. 



B. Symmetry conditions for annihilation 

The dominant electron-positron annihilation processes 
are the 27 and 37 processes. The 27 annihilation rate 
for bound systems is proportional to the probability of 
finding an electron and a positron at the same position 
in a spin-singlet state according to 

T = 47rr^<*|£c^a(r i -r p )|*> 

i 

= 2.018788 x 10 11 ^(<5(r* - r p )) s , (1) 

i 

[H, [25|, [2(| , where the sum is over the electron coordi- 
nates, the 6- function expectation is evaluated in a 3 ,, and 
r is given numerically in s — 1 . The operator Of p is a spin 
projection operator to select spin-singlet states for the i,p 
electron-positron pairs in the wave-function 'J which is 
anti-symmetrized in the electron coordinates. The rate 
constant for the Ps ground state is about 8 x 10 9 s^ 1 . 
Equation ([1]) involves a contact interaction which means 
that the relative angular momentum of the annihilating 
pair (L re i) must be zero [24|. However, electron-positron 
annihilation is possible even for when the relative an- 
gular momentum of the annihilating pair is greater than 
zero. For example the Ps(2p) levels can undergo 27 anni- 
hilation at rates proportional to a 5 and a 6 respectively 
[13, HI]. The rates for the different Ps(2p) levels have 



III. CALCULATION METHODS 

A. The configuration interaction method 

A majority of the calculations in the present pa- 
per wereperformed with a configuration interaction ap- 
proach [23, H3, HH . The CI basis was constructed by let- 
ting the two electrons (particles 1 and 2) and the positron 
(particle 0) form all the possible total angular momen- 
tum Lt — configurations, with the two electrons in a 
spin-triplet state, subject to the selection rules, 

max(4,^ 2 ) < J, (2) 
mm(h,l 2 ) < L- mt , (3) 
(_l)(*>-Mi+40 = _i _ (4) 

In these rules £q, l\ and 1% are respectively the orbital 
angular momenta of the positron and the two electrons. 
We define (E)j to be the energy of the calculation with 
a maximum orbital angular momentum of J. The sin- 
gle particle orbitals were Laguerre Type Orbitals (LTOs) 
with a common exponent chosen for all the orbitals of 
a common i [1^, [30, [3l| . The orbital basis sets for the 
positron and electrons were identical. 

A major technical problem afflicting CI calculations of 
positron-atom interactions is the slow convergence of the 
energy with J [llj Hl|. The J — > 00 energy, (-E)oo, is 
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determined by the use of an asymptotic analysis. The 
successive increments, AEj = (E)j — (E)j_i, to the 
energy can be written as an inverse power series [HI, HI, 
iil, viz 



AEj 



B f 



C f 



(J+I)6 (J+|)7 (J + i)8 (J+ i )8 



(5) 

The first term in the series starts with a power of 6 since 
all the possible couplings of any two of the particles result 
in unnatural parity states [36j . 

The J — > oo limit, has been determined by fitting sets 
of (E)j values to asymptotic series with either 1, 2, 3 or 

4 terms. The coefficients, Ae, Be, Ce and De for the 
4-term expansion are determined at a particular J from 

5 successive energies ((E)j^ 4 , (E)j_ 3 , {E)j- 2 , (E)j-i 
and (E)j). Once the coefficients have been determined it 
is easy to sum the series to oo and obtain the variational 
limit. Application of asymptotic series analysis to helium 
has resulted in CI calculations reproducing the grou nd 
state energy to an accuracy of «10~ 8 hartree [35L l37t| . 

The treatment of the alkalis Li, Na and K requires 
the use of a frozen core approximation. The details of 
this approximation have been discussed in great detail 
elsewhere [2{| |3(| HU , so only the briefest description is 
given here. The model Hamiltonian is initially based on 
a hartree-Fock (HF) wave function for the neutral atom 
ground state. The core orbitals are then frozen. The di- 
rect part of the core potential is attractive for electrons 
and repulsive for the positron. The impact of the direct 
and exchange part of the HF core interactions on the ac- 
tive particles are computed without approximation. One- 
and two-body semi-empirical polarization potentials are 
then added to the potential. The adjustable parameters 
of the core-polarization potential are defined by reference 
to the spectrum of neutral atom [3(| HH . 



B. The stochastic variational method (SVM) 

In the stochastic variational approach an explicitly cor- 
related gaussian (ECG) is constructed by placing the 
particles (electrons and positrons) into Gaussian single 
particle orbitals 

rlY lm {h)e^{-/3^} = JMr.^expl-Arf} , (6) 

and using a 

exp{-aij(ri -r^) 2 } , (7) 

Gaussian correlation function between the ith and jth 
particles. The iV-particle trial function is then 

*Ls(r) - A{[[[y h y h ] li2 y h ] . . ] LMl X sm s 

N 

x n ex p{-& r * 2 }n ex p{-M r .? - r >) 2 }}(8) 



where A is an antisymmetrizer and xsm s 1S t ne spin func- 
tion of the particles. The nonlinear variational parame- 
ters /3i and ctij are selected by an iterative trial and error 
procedure. Full details are given in Ref. [1^, lioj . The 
orbital angular momentum quantum numbers U are re- 
stricted to occupy the lowest possible values consistent 
with the overall symmetry of the state. The spherical 
part of the ECG basis functions effectively allows inter- 
nal angular momentum to be distributed between the 
different parts of the systems. Accordingly, eq. ([8|) im- 
plicitly includes all possible internal symmetries that can 
make a contribution to the energy. This has been verified 
with test calculations. 



IV. RESULTS OF CALCULATIONS 

A. The A~(np 2 3 P C ) states 



TABLE I: The energies of various parent states relevant to 
the structure and energy threshold for the APs unnatural par- 
ity states. The polarizability only allows for np — » kd exci- 
tations since np — > ks excitations cannot occur in the P c 
channel. The energy for the H~(np 2 3 P C ) state was taken 
from Bylicki and Bednarz [gij, while those for the alkali sys- 
tems were from the present CI calculations. 





A(np) 


A(np 2 3 P e ) 


A(np) + Ps(2p) 


Old (do) 


H 


-0.12500 


-0.12535545 


-0.1875 


173.3 


Li 


-0.13023850 


Unbound 


-0.19273850 


142.7 


Na 


-0.11156287 


-0.11382478 


-0.17406287 


302.0 


K 


-0.10018265 


-0.10450418 


-0.16268271 


557.6 



TableUgives the energies of the various parent states of 
the APs systems. These are relevant to the determination 
of the energy thresholds. The energies for the A~(np 2 
3 P C ) states were taken from CI calculations which used 
an exact subset of the basis used for the calculations upon 
the APs system. The energy for the H~(np 2 3 P e ) state 
was taken from a large CI-Hylleraas calculation :23] that 
was converged to eleven significant digits. The results of 
an SVM calculation of this state are detailed in Table HT1 
The dipole polarizabilities listed in Table Q] show an 
interesting correlation between the polarizability of the 
A(np) state and the electron affinity in the A~(np 2 3 P e ) 
channel. The larger the polarizability, the larger the 
binding energy. The Li(2p) level has the smallest polar- 
izability, and is the one atom that is unable to support a 
negative ion in the 3 P e channel. 

This behaviour is reminiscent of the electron affinity 
systematics of the alkaline-earths in the 2 P° channel. The 
ground states of Be and Mg do not have an electron affin- 
ity while those of Ca, Sr and Ba have electron affinities 
that become larger as the atom, and its polarizability, be- 
come larger 41] . The critical polarizability for the alkali 



sequence is somewhere between 142 and 173 a 3 ,. The po- 



larizability of calcium, which just binds an electron with 
an electron affinityof « 7 x 1(T 4 hartree [H [H, Si| , is 



about 160 dg [3 



TABLE II: Properties of the ET( 3 P C ), PsH( 2 ' 4 S°) and 
LiPs( 2,4 S°) ground states. Data are given assuming an in- 
finite nuclear mass (n). All quantities are given in atomic 
units. The magnitude of the binding energy against dissocia- 
tion into the lowest energy fragmentation channels is given by 
e while T+ and T_ represent the positron and electron kinetic 
energy operators. 



changing by about 1%. 



Property 


PsH 


H~ 


LiPs 


N 


400 


400 


1000 


(V) /(T) + 2 


5.1 x 10~ 8 


4.4 x 10~ 9 


5.6 x 10 -5 


E 


-0.188317 


-0.12535545 


-7.472871 


e 


0.000817 


0.00035545 


0.000215 


(T_) 


0.156579554 


0.125355451 




(T+) 


0.031737544 






\ r ne- I 


O.OU ( 


1 1 fiWfil Q 


A "il 77 


(r ne +) 


14.243 




12.991 


(r E -e-) 


14.001 


19.58289 


7.6429 


v"e+ e -) 


12.722 




12.531 


(Vr„.-> 


0.174397 


0.160521 


1.4358 


(V r ne+) 


0.084716 




0.089443 


<l/r.- e -> 


0.089023 


0.0700331 


0.40130 


(l/r e+e -> 


0.101789 




0.096298 


(rL-) 


121.185 


271.2046 


45.959 




247.910 




192.691 


(rl- £ -) 


245.959 


556.893 


92.533 




202.204 




182.800 




B. 


PsH 





1. The configuration interaction method 

The Hamiltonian was diagonalized in a basis con- 
structed from a large number of single particle orbitals, 
including orbitals up to I = 10. There were 20 radial 
basis functions for each I. Note, the symmetry of the 
state prevented the electrons or positrons from occupy- 
ing £ = orbitals. The largest calculation was performed 
with J = 10 and L- lnt = 3 and gave a CI basis dimen- 
sion of 369200. The parameter L- lnt does not have to 
be particularly large since it is mainly concerned with 
electron-electron correlations [3(|. The resulting Hamil- 
tonian matrix was diagonalized with the Davidson algo- 
rithm 46| , and a total of 300 iterations were required for 
the largest calculation. The present calculation is very 
slightly different from that reported in [l2j]. One of the 
I = 10 Laguerre functions in [12] was input with the 
incorrect n. The inclusion of the correct Laguerre func- 
tion resulted in the final binding energy reported in [l2j 




PsH 2 S° 



FIG. 1: The binding energy, e = -((E) + 0.1875), of the 
2 ' 4 S° state of PsH as a function of J. The directly calculated 
energy is shown as the solid line while the J — > oo limits using 
eq. ((5} with 1, 2, 3 or 4 terms are shown as the dashed lines. 
The binding energy of the SVM wave function is also shown. 
The H(2p) + Ps(2p) dissociation threshold is shown as the 
horizontal solid line. . 

The energy of the PsH 2,4 S° state as a function of J 
is given in Table IIIII The calculations only give an en- 
ergy lower than the H(2p) + Ps(2p) threshold of -0.1875 
hartree for J > 9. Figure [1] shows the estimates of (E)^ 
as a function of J. A quick visual examination suggests 
that the extrapolations are converging to a common en- 
ergy which attests to the reliability of the extrapolations 
in J. The impact of the extrapolations is significant since 
they more than double the binding energy. The best CI 
estimate of the binding energy is the four-term extrapo- 
lation at J = 10 listed in Table Hill namely 7.10 x 10 -4 
hartree. The main area where improvement could be 
made is in the dimension of the radial basis. A precur- 
sor to the present CI calculation with 15 LTOs gave an 
extrapolated binding energy of 6.06 x 10 -4 hartree. 

The extrapolations of the other expectation values in 
Table |Hl] were done using eq. ((5|). It should be noted that 
there is no formal justification for the use of eq. (JSJ) for 
expectation values other than the energy, so there is an 
additional degree of uncertainty for these extrapolations. 
In practice, this extra uncertainty is not that significant 
since the finite dimension of the radial basis represents a 
larger source of error. 



2. The stochastic variational method 

For the PsH calculation, the two electrons have been 
placed in l\ = I2 = 1 orbitals and coupled to an L = 1 
state with a total spin of S = 1. The positron is then 
placed into an Iq — 1 orbital and the whole composite 
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TABLE III: The energy of the 2 ' 4 S° state of PsH as a function of J and with L int = 3. The threshold for binding is -0.1875 
hartree. The column n gives the total number of occupied electron orbitals (the number of positron orbitals was the same) 
while Nci gives the total number of configurations. The radial expectation values for the electron, (r e ), and the positron, (r p ) 
are given in ao. The results of the J — > oo extrapolations using eq. (f5| at J = 10 are given. 



J 






n 


Nci 


(E)j 


£ 




(r P ) 


1 






20 


4200 


-0.16755818 


— 0.01994182 


7.08076 


13.10807 


2 






40 


16400 


-0.17938458 


— 0.00811542 


6.95155 


11.88797 


3 






60 


45000 


-0.18327391 


-0.00422609 


7.08870 


11.66425 


4 






80 


85000 


-0.18510516 


-0.00239484 


7.23884 


11.70642 


5 






100 


129200 


-0.18612684 


-0.00137316 


7.37821 


11.82951 


6 






120 


177200 


-0.18675237 


-0.00074763 


7.50464 


11.97767 


7 






140 


225200 


-0.18715897 


-0.00034103 


7 61882 


12.13043 


8 






160 


273200 


-0.18743569 


-0.00006431 


7.72208 


12.27939 


9 






180 


321200 


-0.18763074 


0.00013074 


7.81562 


12.42101 


10 






200 


369200 


-0.18777213 


0.00027213 


7.90047 


12.55387 












J — > oo extrapolations 






1-term 


eq. 








-0.18800504 


0.00050504 


8.04024 


12.77272 


2-term 


eq. 


© 






-0.18811689 


0.00061689 


8.14930 


12.95211 


3-term eq. © 






-0.18817637 


0.00067637 


8.23457 


13.09609 


4-term 


eq. 


© 






-0.18821031 


0.00071031 


8.30116 


13.21039 



is coupled to L = 0. The largest calculation had a to- 
tal of 400 ECGs. A summary of the energy and other 
expectation values is given in Table |TTJ 

The energy of the best SVM wave function was 
-0.188317, yielding a binding energy of 8.17 x 10~ 4 
hartree. The deviation of the wave function from the 
exact virial theorem expectation, ((V)/(T) + 2), was 
5.1 x 10~ 8 . A many body system interacting by purely 
coulombic interactions is known to satisfy ((V) / (T) = 
-2) The SVM binding energy is just over 10% 

larger than the CI energy and should be closer to the 
variational limit. 

One interesting aspect of the 4 S° state is that it is more 
tightly bound than its H(2p 2 ) parent. Its binding energy 
is more that twice as large as the H(2p 2 ) binding energy 
of 3.55 x 10~ 4 hartree j23|. It is also more compact. The 
mean electron distance from the nucleus of (r e ) = 8.86 
ao is smaller than that for H(2p 2 ), namely (r e ) = 11.66 
ao. In short, the addition of the positron has resulted 
in a complex that has a larger binding energy than its 
3-body parent. 

The SVM radial expectation for the positron, (r p ) was 
14.24 a , somewhat larger than the extrapolated CI value 
of 13.67 ao- In the CI calculation the positron is localized 
closer to the nucleus even though the CI wave function is 
less tightly bound. This is a purely computation limita- 
tion, due to the nature of LTO basis which is relatively 
compact. Improving the radial expectation for the CI 
wave function would require an increase in the number 
of radial functions per £. 

The inter-particle correlation function, C(r), is de- 
fined as the probability of finding any pair of parti- 
cles a certain distance apart. The correlation func- 



tions shown in Figures [2] and Figures [3] are consistent 
with a structure consisting of a Ps(2p) complex weakly 
bound to the H(2p) state. Consider an idealized struc- 
ture consisting of a product wave function of the form 
# = $(Ps[2p])$(H[2p])$ Ps ( J R) where $ Ps (i?) is the wave 
function describing the motion of the Ps(2p) center of 
mass. The (p, e~) and (p, e + ) correlation functions aris- 
ing from the Ps(2p) cluster should be the same. There- 
fore, adding the (p, e + ) correlation function to the (p, e~) 
C(r) of H(2p) state should give a correlation function 
that is the same as the actual (p, e~) correlation Figure 
[2] shows a strong degree of resemblance between the ac- 
tual (p, e~) correlation function and that obtained from a 
$(Ps[2p])$(H[2p])$ Ps (i?). Similarly, adding the {e~,e~) 
correlation function to the (e + ,e~) C(r) of Ps(2p) state 
should give a correlation function that is the same as the 
actual (e _ , e~) correlation function. Once again, the two 
curves shown in Figure [3] a degree of similarity. 

The energies of the finite mass variants of PsH have 
also been determined. The energies of Ps 1 !!, Ps 2 H, 
and Ps 3 H arc -0.1882398, -0.1882784 and -0.1882913 
hartree respectively. The binding energies are 8.078 x 
10~ 4 , 8.124 x 10~ 4 , 8.140 x 10~ 4 , hartree respectively. 



C. LiPs 

The 2 S° state of LiPs is a very unusual state in that it 
is a Borromean state [HI, [3 0, EH] ■ This is because all 
the possible 3-body parent states, namely the 3 P° states 
of Li - , e + Li or Ps~, are themselves unstable. 

The Ps~ ion has been thoroughly investigated and does 
not possess a stable 3 P e state [471. |48| . 
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FIG. 2: The correlation functions for the (p,e~) and (p, e + ) 
particles of PsH. Also shown is a correlation function obtained 
by adding the (p, e + ) correlation function to the (p,e~) cor- 
relation function of the H(2p) state. 



0.20 




whether the e + H state is stable in the 3 P e channel, and 
once again there was no indication of a bound state. 

The calculations upon LiPs were very similar in scope 
and scale to those carried out upon PsH although the 
calculations were taken to J = 11 in order to have an 
explicit calculation that gave binding. The sequence of 
CI energies and other expectation values as a function of 
J are given in Table ITVl The binding energy ej is defined 
as sj = ~({E) + 0.19273850). 

Figure[3]depicts the binding energy and extrapolations 
as a function of J. Only for the J = 11 basis has (e)j 
crossed the threshold for binding. 

The most reliable estimates of the energy is that given 
after the 4-term extrapolation is used to determine the 
J — > oo limit of the binding energy. The different curves 
in Figure [4] tend to be closer together as the number of 
terms in the extrapolation increase. The binding energy 
of 4.42 x 10~ 4 hartree is just over half that of the PsH 
state. 



1.0 



S 



0.0 



-1.0 < L 



LiPs 2 S° 



£„, 3 -term 
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FIG. 3: The correlation functions for the (e~,e~) and 
(e~,e + ) particles of PsH. Also shown is a correlation func- 
tion obtained by adding the (e~,e~) correlation function to 
the (e~,e + ) correlation function of the Ps(2p) state. 



FIG. 4: The binding energy of the 2,4 S° state of LiPs as a 
function of J. The directly calculated binding energy is shown 
as the solid line while the J — > oo limits using eq. ([5| are 
shown as the dashed lines. The Li(2p) + Ps(2p) dissociation 
threshold is shown as the horizontal solid line, lines. 



Similarly, the Li~ ion is believed not to have a stable 
3 P C state [i^, We have also performed some very 

large CI calculations upon the Li - ion and these calcu- 
lations gave no indication of a bound state in the 3 P e 
symmetry. 

Finally, the e + Li system is also not stable in the 3 P e 
channel. Once again a very large CI calculation has been 
performed and once again there was no indication of a 
bound state. Further, some calculations of e + -Li scat- 
tering in the 3 P e channel also gave no sign of a bound 
state. The polarizabilities given in Table U also indicate 
that it should be easier to bind a positron to the H(2p) 
state than the Li(2p). The SVM was also used to check 



The positron can annihilate with the core electrons 
via the 27 process since the symmetry considerations are 
irrelevant here. However, the annihilation rate of r coro w 
4 x 10 4 s _1 is small because the positron cannot occupy 
a £ = orbital. 

The mean positron-nucleus distance of (r p ) = 13.5 do 
for the CI wave function was almost the same as the 
CI wave function estimate for PsH despite the smaller 
binding energy. Part of the reason for this lies in the 
LTO basis sets which were almost identical for the two 
atoms. The finite range of the LTO basis could be acting 
to artificially confine the positron. However, it must be 
remembered that the asymptotic Ps(2p) cluster will also 
be confined by the L — 1 centrifugal barrier. 
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TABLE IV: Results of CI calculations for the 1 S° state of LiPs for a series of J, with fixed L- m t = 3. The 3-body energy of the 
system, relative to the energy of the Li + core, is denoted by E (in hartree). The threshold for binding is -0.19273850 hartree, 
and e gives the binding energy (in hartree) against dissociation into Ps(2p) + Li(2p). The core annihilation rate in units of 
s _1 is given in the Y c column. The numbers in square brackets indicate powers of 10. Other aspects of the table design are 
identical to those of Mil 



J 


n 


N C i 


E 




(r r ) 


(r„) 


r c 


1 


20 


4200 


—0 17291946 


—0 01981904 


6 86606 


12.84274 


4.5441 [4] 


2 


40 


16400 


—0 18449498 


—0 00824351 


6 75226 


11 69309 


6 5727F41 


3 


60 


45000 


-0.18828973 


-0.00444877 


6.89326 


11.50143 


6.5611[4] 


4 


80 


85000 


-0.19007657 


-0.00266192 


7.04675 


11.56062 


6.1767[4] 


5 


100 


129200 


-0.19107655 


-0.00166195 


7.18989 


11.69580 


5.7955[4] 


6 


120 


177200 


-0.19169133 


-0.00104717 


7.32070 


11.85451 


5.4761 [4] 


7 


140 


225200 


-0.19209300 


-0.00064550 


7.43986 


12.01755 


5.2173[4] 


8 


160 


273200 


-0.19236789 


-0.00037061 


7.54860 


12.17700 


5. 0080 [4] 


9 


180 


321200 


-0.19256282 


-0.00017568 


7.64806 


12.32954 


4.8376[4] 


10 


200 


369200 


-0.19270505 


-0.00003345 


7.73914 


12.47367 


4.6972 [4] 


11 


220 


417200 


-0.19281127 


0.00007278 


7.82273 


12.60906 


4.5806[4] 








J -> 


oo extrapolations 








1-term eq. 


© 




-0.19300706 


0.000268567 


7.97681 


12.85861 


4.3656[4] 


2-term eq. 


© 




-0.19310120 


0.000362688 


8.09858 


13.06361 


4.2226[4] 


3-term eq. 


© 




-0.19315143 


0.000412934 


8.19533 


13.23002 


4. 1331 [4] 


4-term eq. 


© 




-0.19318032 


0.000441825 


8.27254 


13.36466 


4.0698[4] 



1. The stochastic variational method 

For the SVM, in the LiPs case, the first two electrons 
are placed in the l\ = li — orbits and their spins 
are coupled to zero, the next two electrons are in the 
I3 =£4 = 1 orbits (with the total angular momentum 
coupled to 1) and their spins are coupled to 1. Finally, 
the positron is placed in an Iq = 1 orbit and the total 
orbital angular momentum is coupled to 1. 

The threshold for binding is the Li(2p) (E = 
—7.4101565 hartree [5l|) plus the Ps(2p) energy. So the 
energy threshold for an absolute variational proof of bind- 
ing is at —7.4726565 hartree. The energy and expectation 
values of the SVM LiPs wave functions are listed in Table 
ILT1 The best variational energy was —7.472871 hartree, 
equivalent to a binding energy of 2.15 x 10~ 4 hartree. 
The energy optimization was not fully completed and the 
binding energy of the CI calculation is probably more re- 
liable. The primary purpose of the SVM calculation was 
to give an absolute proof that the unnatural parity state 
of LiPs was electronically stable. 



D. NaPs 



the model potential were -0.18885491 and -0.11156287 
hartree. The experimental binding energies are 
-0.188858 and -0.111547 hartree respectively 
Electronic stability requires a total 3-body energy of 
—0.17406287 hartree and the binding energy ej is de- 
fined as ej = -((E) + 0.17406287). The energy of the 
3 P e excited state of Na" is -0.11342529 hartree, i.e the 
Na(3p) has an electron affinity of 0.002262 hartree with 
respect to attaching an electron to the 3 P e state. This is 
reasonably close to the original value of Norcross, 0.00228 
hartree [49| |. 

Table [V] gives the energies and radial expectation val- 
ues as a function of J while figure [5] shows the variation 
of £oo as a function of J. Once again the three and four 
term extrapolations seem to be converging to a common 
energy. In this case the J — > 00 correction increases the 
binding energy by about 40 % from 10.87 x 10~ 4 hartree 
to 15.14 x 10 -4 hartree. The binding energy of the NaPs 
unnatural parity state is about twice as large as that of 
PsH. 

The annihilation rate with the core electrons was 
r core w 1.5 x 10 5 s _1 is small. Although this is 3 times 
larger than r core for LiPs, in absolute terms the annihi- 
lation rate is still small. 



The calculations upon NaPs were very similar in scope 
and scale to those carried out upon LiPs. About the only 
difference was that an extra I — 1 orbital was added to 
the electron basis. 

The energies of the Na(3s) and Na(3p) states in 



E. KPs 

The calculations upon KPs were very similar in scope 
and scale to those carried out upon LiPs. About the only 
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TABLE V: The energy of the 2,4 S° state of NaPs as a function of J. Energies are given relative to that of the Na + core while 
the threshold for binding is —0.17406287 hartree. The column n_ gives the total number of occupied electron orbitals, while 
n+ gives the number of positron orbitals. Other aspects of the table design are identical to those of II VI 



J n_ 


n + 


± »CJ 






\'e/ 


\' pi 


r 

r c 


1 91 


9D 




fl 1 5^785fiQ 
— u.ioo i oouy 


n 0909771 8 


7 74987 


1 "3 SSfiflS 

lO.ODUUO 


1 .4441 1 O ] 


9 41 




1 799f) 


fl 1 fifil 4955 


— u.uu j yiUOi 


J . OOO lO 


1 9 4RQ5D 


9 9' : i' : !Q[5l 


3 fil 


Kfi 
uu 


4fi99n 


D 1 70' : !^951 


n nn ,: !7^n^5 

u.uuo / ouoo 


7 71 D^Q 


1 9 1 RR54 


9 9Q48f5l 


4 81 


80 


86620 


-0.17231600 


-0.00174687 


7.85240 


12.16161 


2.1849[5] 


5 101 


100 


131220 


—0 1 7341 763 


—0 00064523 


7 98365 


1 2 25289 


2 0625 [51 


6 121 


120 


1 79620 


—0 1 7408661 


00002375 


8 1 0085 


1 2 37548 


1 9577f5l 


7 141 


140 


228020 


—0 17451658 


00045371 


8 20440 

( ) . Z- V 7 i tv' 


1 2 50450 


1 8733 [51 


8 161 


160 


276420 


—0 1 7480552 


00074266 


8 29583 


12 62983 


1 8060 [51 


9 181 


180 


324820 


-0.17500636 


0.00094349 


8.37645 


12.74713 


1.7523 [5] 


10 201 


200 


373220 


-0.17514972 


0.00108685 


8.44734 


12.85449 


1.7093 [5] 








J — » (X) 


extrapolations 








1-term eq. ([5]) 






-0.17538587 


0.00132300 


8.56412 


13.03133 


1.63839[5] 


2-term eq. (J5]) 






-0.17549381 


0.00143094 


8.65082 


13.17118 


1.59424[5] 


3-term eq. ([5]) 






-0.17554787 


0.00148500 


8.71500 


13.27807 


1.56638[5] 


4-term eq. §5§ 






-0.17557663 


0.00151376 


8.76225 


13.35833 


1.54841 [5] 



TABLE VI: The energy of the 2 ' 4 S° states of KPs as a function of J. The threshold for binding is -0.16268265 hartree and 
and the energies are given relative to that of the K + core. Other aspects of the table design are identical to those of Table [V] 



J 71- 


n+ 


Nat 


E 


E 


(r.) 


(r P ) 


r c 


1 21 


20 


4620 


-0.14168301 


-0.02099965 


8.49503 


14.71924 


3.4488[5] 


2 41 


40 


17220 


-0.15457885 


-0.00810380 


8.29108 


13.13410 


5.4508[5] 


3 61 


60 


46220 


-0.15914469 


-0.00353796 


8.39896 


12.70578 


5.7120[5] 


4 81 


80 


86620 


-0.16132591 


-0.00135675 


8.53370 


12.64979 


5.4824[5] 


5 101 


100 


131220 


-0.16253700 


-0.00014566 


8.66023 


12.71399 


5.1917[5] 


6 121 


120 


179620 


-0.16326995 


0.00058729 


8.77291 


12.81839 


4.9358[5] 


7 141 


140 


228020 


-0.16373804 


0.00105538 


8.87154 


12.93315 


4.7280[5] 


8 161 


160 


276420 


-0.16405038 


0.00136772 


8.95749 


13.04566 


4.5629[5] 


9 181 


180 


324820 


-0.16426580 


0.00158314 


9.03227 


13.15085 


4.4319[5] 


10 201 


200 


373220 


-0.16441830 


0.00173564 


9.09676 


13.24589 


4.3282[5] 








J — > CO 


extrapolations 








1-term eq. ([5]) 






-0.16466950 


0.00198684 


9.20300 


13.40245 


4.1575[5] 


2-term eq. ([5]) 






-0.16478114 


0.00209848 


9.27933 


13.52387 


4.0541 [5] 


3-term eq. © 






-0.16483516 


0.00215250 


9.33339 


13.61318 


3.9924[5] 


4-term eq. (J5]) 






-0.16486273 


0.00218008 


9.37076 


13.67636 


3.9562[5] 



difference was that an extra £ — 1 orbital was added to 
the electron basis. 

The energies of the K(4s) and K(4p) states in 
the model potential were -0.159520 and -0.10018265 
hartree. The experimental binding energies are 
-0.159516 and -0.100176 hartree respectively [13]. 
Electronic stability requires a total 3-body energy of 
—0.16268265 hartree and the binding energy sj is defined 
as Ej = -((E) + 0.16268265). The energy of the 3 P C ex- 
cited state of K" is -0.104498 hartree, i.e the K(4p) has 
an electron affinity of 0.004322 hartree with respect to 



attaching an electron to the 3 P e state. This is close to 
the original value of Norcross, 0.00437 hartree [49| . 

Table [VTI gives the energies and radial expectation val- 
ues as a function of J while figure[S]shows the variation of 
£oo as a function of J. The three and four term extrapola- 
tions seem to be converging to a common energy. In this 
case the J — > oo corrections increase the binding energy 
by about 20% from 17.36 x 10" 4 hartree to 21.80 x 10~ 4 
hartree. The KPs system has the largest binding energy 
of all the systems considered in this paper. 
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E 



e„ 4-term 




NaPs 2 S° 



10 J 



FIG. 5: The binding energy of the 2 ' 4 S° state of NaPs as a 
function of J. The directly calculated binding energy is shown 
as the solid line while the J — > oo limits using eq. ([5]) with 
1, 2 or 3 terms are shown as the dashed lines. The Na(3p) + 
Ps(2p) dissociation threshold is shown as the horizontal solid 
line. 
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FIG. 6: The binding energy of the 2,4 S° state of KPs as a 
function of J. The directly calculated binding energy is shown 
as the solid line while the J — > oo limits using eq. ([5]) are 
shown as the dashed lines. The K(4p) + Ps(2p) dissociation 
threshold is shown as the horizontal solid line, lines. 



unusual of the systems is LiPs since the 2 ' 4 S° states are 
of Borromean type. The sequence of calculations suggest 
that there would also exist unnatural parity 2,4 S° com- 
plexes of RbPs and CsPs; and most likely they would 
have binding energies larger than KPs. 

Due to their low binding energies, these systems can 
be expected to have a structure composed of an Ps(2p) 
cluster loosely bound to an atomic X(np) excited state. 
This has been confirmed by the correlation functions for 
PsH which were computed using the SVM. 

Although these complexes are electronically stable and 
decay very slowly by electron-positron annihilation there 
are other decay processes that act to shorten the lifetime. 
These complexes can emit a photon, decaying to a state 
of 2 - 4 P c symmetry. For example, a Ps(np) fragment in 
the complex can emit a photon decaying to a Ps(ls) type 
fragment. The Ps(ls) fragment could then annihilate by 
the 27 or 37 process. In addition, the resulting 2 > 4 P C 
state could also decay by auto-ionization. The lifetime 
of these states can be expected to be comparable to the 
lifetime of the fragments against single photon decay, e.g. 
X(2p) — -> X(ls). So the overall lifetimes of the states can 
be expected to be of order 10~ 8 - 10~ 9 seconds. 

It is unlikely that any of these complexes will be iden- 
tified in the laboratory in the near future. The formation 
of positronic compounds is known to be notoriously dif- 
ficult [53j]. That these states are unnatural parity states 
compounds the difficulty since such states are not read- 
ily formed in normal collision systems. For example, the 

3 P° ion states [H SHHHI that C ° Uld S6rVe aS suitable 
parents have never been identified in the laboratory. 

Besides the PsH and APs systems, there are other re- 
lated physical systems that could have unnatural parity 
bound states. For example, there is the possible exis- 
tence of a new bi-exciton excited state [5J|. While the 
Ps~ ion might not have a stable 3 P e state, it is known 
that the (M+, e~, e~)ion is stable for M + /m e < 0.4047 
and M + /m e > 16.8 0, El]. It could be expected that 
a bi-exciton state, (e~,e~,/i, h), with 1 ' 3 > 5 S° symmetry 
would be electronically stable when the mass ratios make 
the 3 P e state of the charged exciton (e~~,e~,h) stable. 
The system might also exhibit Borromean binding, there 
might be a bound bi-exciton state even though neither 
of the 3 P C (e~ ,e~ ,h) or (e~,h,h) states was stable. 



SUMMARY 



Acknowledgments 



A number of PsX systems (X = H, Li, Na and K) are 
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